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Abstract 

We discuss a recent method proposed by Kryukov and Walton to address boundary- value 
problems in the context of deformation quantization. We compare their method with our own 
approach and establish a connection between the two formalisms. 
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^In a recent paper [1] Kryukov and Walton have addressed the problem of infinite walls in the context 
Q©f deformation quantization (see also [2]). They proposed a method alternative to ours [3, 4]. In 
j>!this comment, we would like to compare the two approaches and prove their equivalence under some 
'^simplifying conditions. 

Without loss of generality, we shall henceforth consider a one-dimensional system on the half-line 
< 0, whose state is described by a wavefunction satisfying the Dirichlet boundary condition 

V(o) = o. (i) 

As in [1] we shall take units such that % = 2m = 1. For simplicity, we shall restrict our analysis to 
the free particle. 

In standard quantum mechanics, one solves the time independent Schrodinger equation, 

-ip"(x) = Ei/j(x), (2) 
on the whole line, imposes the condition (1) and confines the wavefunction thereafter: 

(f)(x) = 9(-x)ip(x). (3) 
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In our work we proved that the corresponding Wigner function 

1 r+oc 1 r—x 

f{x,p) = - dy e- 2tpy <P*(x - y)4>{x + y) = 8(-x)- / dy e~ 2ipy ^*{x - y)^{x + y) (4) 

does not obey the expected *-genvalue equation: 

p 2 *f(x,p) = Ef(x,p). (5) 
Since this equation is the Weyl transform of the eigenvalue equation: 

p 2 |0x0|=£|0><0| (6) 

one is led to the conclusion that either the Weyl transform or the eigenvalue equation (6) fail to 
be valid in the confined case. Let us discuss the first possibility in some detail. Standard operator 
quantum mechanics provides two main (non-equivalent) formulations of the confined particle subject 
to Dirichlet boundary conditions, [5, 6]. In the first one the confined system is formulated on 
the Hilbert space L 2 (I), where I is the support of the wave function (typically the half-line or a 
compact interval). In this case the observables' algebra is no longer the Heisenberg algebra (notice 
that there is no self-adjoint representation of the Heisenberg algebra on L 2 (I) which is compatible 
with Dirichlet boundary conditions) [6]. Since the Weyl transform only applies to the Heisenberg 
algebra, it cannot be used in this case. Another possibility is to formulate the system on the Hilbert 
space L 2 (1R) (i.e. acknowledging the rest of the real line) in which case the observable's algebra is 
the standard Heisenberg algebra and the confinement of the particle is a consequence of a suitable 
potential added to the Hamiltonian. In this case we say that the confinement is dynamical and the 
Weyl map is valid. The problem with this approach is that there is no regular potential able to 
confine completely the particle. We then have two possibilities. We either i) provide a family of 
regular potentials displaying the confined solution in the limit case or ii) we introduce a non-regular, 
distributional potential. 

In both papers [1, 3] the authors took the dynamical confining point of view. The two approaches 
differ precisely in the way the boundary potential is modelled. In ref. [1] the authors consider a 
family of regular exponential potentials V a (x) = e 2ax and were left with the energy eigenvalue 
equation: 

(p 2 + V a (x)) Ux) = E U%) (7) 

whose solutions satisfy: lim a ^. 00 £ a (x) = <p(x). 
Applying the Weyl map to (7) we get: 



(p 2 + e 2ax ) -k f a (x,p) = Ef a (x,p) 
f a (x,p) * (p 2 + e 2ax ) = Ef a (x,p) 



(8) 



whose solutions (as pointed out in [1, 2]) also display a well defined limit a — > oo: 

Jim f a (x,p) = f(x,p) (9) 

This means that if W is the Weyl map then lim^oo W(C, a ) = W^liniQ^oo £ a ) a property that 
reinforces our previous argument that the Weyl map properly extends to the confined case provided 
"the rest of the real line" is taken into account. Unfortunately, the pair of equations (8) (as well 
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as eq.(7)) do not make sense in the limit a — > oo given the divergent nature of the potential. 
The authors of [1] were able to manipulate eq.(8) and derive a new equation, alternative (but not 
equivalent) to (8) and displaying the following form in the limit a — > oo: 

= ^d 4 x p(x,p) + l -{p 2 + E)d 2 x p(x,p) + (p 4 - 2Ep 2 + E 2 )p(x,p). (10) 

Moreover they proved that this equation is equivalent to: 

(p 2 - E)* p(x,p)*(p 2 - E) =0 (11) 

and displays a solution which is related to the confined Wigner function (4) by f(x,p) = 9(—x)p(x,p). 

Our purpose in [3] was also to derive a modified *-genvalue equation satisfied by (4). However, 
we followed a completely different procedure: First we proved that <p(x) (cf.(3)) obeys the following 
equation on the whole line: 

-4>"{x) = E<j>(x) - 8'_{x)4>{x), (12) 
where the generalized distribution 5- (x) is defined by: 

J dx S-(x)t(x) = lim + J dx S(x)t(x — e), (13) 

for suitable test distributions t(x) [3, 4]. 

Conversely, it was shown that, if is a solution of (12) and obeys the boundary condition 

0(0-) = 0, (14) 

then 4> takes the form (3), where ip is a solution of (1,2). Eq.(12) defines a modified Hamiltonian: 

H = p 2 + 5'_(x). (15) 

which displays a distributional potential. Applying the Weyl map we get the new *-genvalue 
equations: 

(p 2 + 6'_(x))*f(x,p) = Ef(x,p), f(x,p)*(p 2 + 5'_(x))=Ef(x,p). (16) 

which constitute a system of integral-differential equations for the Wigner function. Notice that 
the integral term comes from the delta distribution: S'_(x) * f(x,p) = lim e ^ + ^ / dkike lkx f(x — 
e,p-k/2),[3}. 

We also derived the boundary conditions to be satisfied by f(x,p). Applying the Weyl map to 
eq.(14) we got: 

lim [dp f(x,p)=0, (17) 

Finally, we proved the converse result in the case of a non-degenerate spectrum. If f(x,p) is a 
solution of (16), and satisfies the boundary condition (17) then f(x,p) takes the form (4) where the 
corresponding wavefunction is a solution of (12,14). 

Several comments are now in order: 
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1) The Kryukov- Walton approach [1] displays several interesting features. It formulates the 
eigenvalue problem in terms of a regular differential equation and it proves that through a suitable 
manipulation the limit a — > oo can be consistently performed in the deformation context. Contrary 
to our own approach one is not required to deal with a differential equation with distributional 
coefficients coming from a distributional potential. This avoids having to deal, in the deforma- 
tion context, with an integral-differential formulation of the eigenvalue problem. Furthermore the 
Kryukov- Walton approach is closer to the spirit of standard local boundary- value problems. Indeed, 
in the latter context one solves a local differential equation, imposes suitable boundary conditions 
and confines the domain thereafter (cf.(2,3)). 

2) On the other hand in [1] one loses contact with the standard operator formulation of the eigen- 
value problem. Strictly speaking eq.(10) cannot be regarded as the Weyl transform of a confined 
version of the eigenvalue equation in standard operator quantum mechanics. In fact its solutions are 
not the Weyl transform of the confined eigenfunctions. Instead, some of its solutions are related to 
the Weyl transform of the confined eigenfunctions by f(x,p) = 9(x)p(x,p). However, equation (10) 
has many solutions. Kryukov and Walton did not provide the conditions that select the ones re- 
lated to the confined stargenfunctions. This is related to the problem of extending Baker's converse 
construction for this formulation (i.e. the proof that if p{x,p) satisfies eq.(10) and some possible 
extra conditions still to be determined then 9(x)p(x,p) is the Weyl transform of the correspond- 
ing confined energy eigenstate). It is also not clear how can eq.(10) be extended to non-diagonal 
stargenfunctions nor if a similar approach can be used to derive the dynamics of the confined system. 

3) In [3] the confined eigenvalue problem is not formulated in terms of a regular potential nor as 
the limit of a family of regular potentials. Instead the boundary potential is a distribution yielding 
a formulation that makes sense quantum mechanically (i.e. it leads to the expected solutions) but 
which is difficult to see how it can be derived by quantizing a classical formulation of the confined 
particle. In this sense the modification of the Hamiltonian introduced in [3] is not derived from 
first principles. Furthermore, the distributional boundary potential added to the Hamiltonian is 
proportional to h 2 contrary to what is standard in the deformation context where typically the 
Hamiltonian is not modified by quantum corrections (i.e. it is just the classical Hamiltonian). 

4) On the other hand in [3] the confined stargenvalue equation is just the Weyl transform of a 
new confined eigenvalue equation valid in standard operator quantum mechanics. Several important 
results follow from this property. In [3] we were able to: i) supply the boundary conditions to be 
satisfied by the confined stargenfunctions (this are just the Weyl transform of the Dirichlet bound- 
ary conditions on the wave function), ii) extend Baker's converse construction to the confined case, 
iii) extend the formalism to non-diagonal stargenfunctions and iv) generate a dynamical (Moyal) 
equation for the confined particle using exactly the same modification of the Hamiltonian. 

To summarize: The method of Kryukov and Walton can be derived from first principles (i.e. by 
quantizing a regular potential) while our own approach is mathematically more complete. In some 
sense the two approaches complement each other. 

Hence an interesting result would be the proof of the equivalence of the two formulations. In 
fact the original aim of the authors of [1] was "to derive the prescription of [3] from first principles". 
However, they found that their "results do not relate easily to the proposal of [3]". 
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In this paper we will prove the equivalence of the two formulations for the simplest case of a 
free particle confined to the negative half-line. We will also provide the necessary boundary and 
kinematical conditions for the approach of [1] (thus proving Baker's converse construction for the 
simplest case mentioned above). Unfortunately, we were unable to extend these results to the time 
dependent case. We discuss this issue in section 3. 

2 Equivalence of the two formulations 

We shall now try to explain how equations (10) and (16) are related. The crux of the matter is the 
fact that (16) are equations for the Wigner function f(x,p) whereas (10) is an equation for p(x,p). 
The two quantities are related by: 

f(x,p) = 9(-x)p(x,p). (18) 

Following the method of [1], we get from eq.(16): 

H*f*H = E 2 f ^p 2 *f*p 2 + p 2 *f*5'_( x ) + 5i(x)*f*p 2 + 5'_(x) *f*6'_(x) = E 2 f. (19) 

From (16), we also have: 

S'_(x) * / = Ef - p 2 * /, / * S'_(x) =Ef-f*p 2 . (20) 

Substituting (20) in (19) we get: 

2ERe(p 2 * /) - p 2 * / * p 2 + 5'_(x) * / * 5'_(x) = E 2 f. (21) 

Now, as stated before, the solution of (16,17) takes the form (4) where the wavefunction obeys a 
Dirichlet boundary condition. In [3] we proved that f(x,p) must then satisfy: 

f(0~,p) = d x f(0',p) = d 2 J(0~,p) = 0, VpeM, (22) 

a result that follows from eq. (4) and which implies: 

p(0, p) = d x p(0, p) = d 2 x p(0, p) = 0, \/p G IR. (23) 

From eq.(18) it follows that: 

d x f(x,p) = -5(x)p(x,p) + 9(-x)d x p(x,p) = 9(-x)d x p(x,p) (24) 

where eq.(23) was taken into account. Likewise: 

' d 2 f(x,p) = 9(-x)d 2 p(x,p) 

< dlf(x,p) = 9(-x)dlp(x,p) (25) 
[ %f(x,p) = 9(-x)dip(x,p) - 5(x)d 3 xP (0,p) 

On the other hand, a simple calculation following [3] leads to: 

S'_(x)*f(x,p)*6'_(x) = ±6(x) |^'(0)| 2 . (26) 



From (18,21,24,25,26), we get: 

(2Ep 2 -p*)6(-x)p - l(E + p*)e(-x)dlp - -L [9(-x)dip(x,p) - 6(x)%p(0,p)] + 

(27) 

+ ±S(x)\^(0)\ 2 = E"9(-x)p(x,p). 
Finally, if we compute d^p(0,p), using (4), we obtain: 

dip{^p) = -lwm 2 . (28) 

Substituting (28) in (27) we obtain precisely (10). This means that our equation (16) together with 
the boundary condition (17) implies the main result of Kryukov and Walton (10). 

We will now argue that the two approaches are in fact equivalent, provided we require that the 
solution of (10) be subjected to the boundary conditions (17,23) and to a pure state condition to 
be stated below [7]. 

We start by looking for solutions of (10) of the form e ir ^ x . From (10) we obtain the 4 solutions: 
flip) — fliv) = 2^G»+ y/~E), r 3 (p) = r%(p) = 2i(p — \fE). Therefore, the most general real solution 
of (10) is: 

p(x,p) = A(p)e 2i ( p+ ^ x + A*(p)e- 2l( < p+ ^ x + B(p)e 2i(p ~^ )x + B*(p) e - 2l ^^ x . (29) 

Now, imposing the boundary conditions (23), which are necessary (albeit not sufficient) to ensure 
that V(0) = [3], we get: 

f sin \2xip + VE)} sin \2x(p - v 7 ^)! ] 

'<'■"> = 1 p+ ve ~ [ p -ve (30) 

where N(p) is an arbitrary real function of p. Notice that if we choose 

Nip) oc -, (31) 
p 

we obtain (up to a multiplicative real constant) the physical solution (eq.(30) of ref.fl]). To derive 
(31) we still have to impose a pure state condition and condition (17). In ref.[7] we proved that if 

E{y,p) = jdx e ix vfix,p), (32) 

obeys the nonlinear partial differential equation 

d 2 1 d 2 

lnE(y,p) = -— -m£(y,p), (33) 



dy 2 " ' 1 Op 

then fix,p) is a pure state Wigner function. Let us apply this condition to fix,p) = 9(—x)p(x,p), 
where pix,p) is given by (30): 

lnE(y,p) =ln[JV(p)p]+lni; 1 ( 2 /,p), (34) 
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where: 



txy 



- / * - \ ^ (35) 



From the comments made after eq.(31), we already know that lnEi(y,p) satisfies (33). We are left 
with: 

— In \pN(p)} = 0, (36) 



the solution of which is 



e 



ap 



N(p) oc — , (37) 



where a is an arbitrary real constant. We then get: 

»+oo 



f(x,p) = 6(-x)p(x,p) oc / dy e 2iyp <p*(x - y)y{x + y), (38) 

J — CO 

where 

<p(x + ia/2) = 9(-x) (e lx ^ - e~ ixV ^) . (39) 

Finally, if we require f(x,p) to satisfy (17), we get a = 0, which is the desired solution (and which 
also avoids the continuation of the Heaviside function into the complex plane). 

So, it seems that the problem can be formulated equivalently in the two following distinct ways. 
Either f(x,p) is a solution of (16,17) or, alternatively, f(x,p) = 9(—x)p(x,p), where p(x,p) is a 
solution of (10,17,23,32,33). The approach of Kryukov and Walton (as our own) remains valid if 
an additional non-singular potential is added to the Hamiltonian. We expect the equivalence of the 
two formulations to be extendable to this more general case. However this is still lacking a proper 
proof. 



3 Discussion 

From the previous analysis it seems that the two methods are equivalent provided suitable boundary 
and kinematical conditions are imposed. At least this seems to be the case for the time-independent 
stargenfunctions of non-degenerate Hamiltonians. By proving this equivalence for the time indepen- 
dent case we also provided a derivation of the approach of [3] from first principles. We also supplied 
the boundary and kinematical conditions selecting the confined stargenfunctions in the method of 
[1] thus extending Baker's converse construction for this approach. Notice that these results still 
have to be extended to the case where an additional regular potential is added. 

On the other hand, the Kryukov- Walton method still lacks a suitable extension to the non- 
diagonal and to the time dependent cases. Unfortunately we think that a time-dependent equation 
is difficult to obtain with this approach. Indeed, the key to deriving eq.(10) is the fact that the 
stargenvalue equation has both a real and an imaginary part, which allow us to express d x p and d^p 
in terms of the quantities p(x,p ± ia) and obtain eq. (37) of ref.[l]. However, the Moyal equation, 

d 

ih—p(x,p,t) = H(x,p)*p(x,p,t) - p(x,p,t)*H(x,p), (40) 
only has a real part. Hence, the previous method cannot be applied straightforwardly. 
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Alternatively, we may follow the method presented here and start from the confined Moyal 
equation (eqs. (41,42) below). It is conceivable that one may eliminate the delta functions appearing 
in this equation. In [3] we proved that for the free particle it reads: 



d t p + 2pd x p = K,(x,p,t), 



(41) 



where 



K{x,p,t) = — e~ 2ip V*(0,t)^(2a;,t) -e 2ip V(2x,t)V'(0,t) 



(42) 



Now suppose that the latter expression obeys the differential equation: 



d t K =V(d x ,d p )lC, 



(43) 



where V is some differential operator. Then, from (41,43) we get: 



(d t -V(d x ,d p ))(d t p + 2pd xP ) = 0. 



(44) 



In any case, even if it were possible to obtain such an equation, the previous analysis reveals that 
this will most likely be a second-(or higher-)order equation with respect to time. This in turn spoils 
one of the most important features of deformation quantization, namely that the time-evolution 
should be a deformation of the classical Liouville equation. 
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